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We consider the Hamiltonian for neutrino oscillations in matter in the case of arbitrary potentials 
including off-diagonal complex terms. We derive the expressions for the corresponding Hamiltonian 
in the basis of the instantaneous eigenstates in matter, in terms of quantities one can derive from 
the flavor basis Hamiltonian and its derivative, for an arbitrary number of neutrino flavors. We 
make our expressions explicit for the two-neutrino flavor case and apply our results to the neutrino 
propagation in core-collapse supernovae where the Hamiltonian includes both coupling to matter and 
to neutrinos. We show that the neutrino flavour evolution depends on the mixing matrix derivatives 
involving not only the derivative of the matter mixing angles but also of the phases. In particular, 
we point out the important role of the phase derivatives, that appear due to the neutrino-neutrino 
interaction, and show how it can cause an oscillating degeneracy between the diagonal elements of 
the Hamiltonian in the basis of the eigenstates in matter. Our results also reveal, that the end of the 
synchronization regime is due to a rapid increase of the phase derivative, and identify the condition 
to be fulfilled for the onset of bipolar oscillations involving both the off-diagonal neutrino-neutrino 
interaction contributions and the vacuum terms. 

PACS numbers: 14.60.Pq, 97.60.Bw 
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I. INTRODUCTION 

Successful core-collapse supernova explosions seem to be 
at reach in the next decade. Currently, sophisticated 
two and three dimensional simulations comprise convec- 
tion, hydrodynamic instabilities (in particular the SASI 
mode), realistic nuclear networks and neutrino transport. 
It appears that the explosion mechanism is being delin- 
eated for a number of progenitor masses [l|, 0] • 

Since the discovery of neutrino oscillations, neutrino 
experiments have determined most of the parameters of 
the Maki-Nakagawa-Sakata-Pontccorvo matrix Q , relat- 
ing the flavor to the mass basis. The focus of next gener- 
ation experiments is to measure the third mixing angle, 
the (Majorana and Dirac) CP violating phases, the ab- 
solute mass scale and hierarchy, the Dirac versus Majo- 
rana neutrino nature The experimental progress has 
numerous implications, e.g. for the neutrino flavor con- 
version in media. It is now clear that the solar neutrino 
deficit is due to Mikhcev-Smirnov-Wolfcnstcin (MSW) 
effect [1, @ - the resonant conversion produced by the 
neutrino interactions while travelling through matter. 

The investigations of the neutrino flavor conversion 
in core-collapse supernovae have revealed an unexpected 
complexity compared to the solar case. This is due to 
the realization of calculations including the neutrino- 
neutrino interaction 0, Q , using matter density profiles 
with shocks (9l-[l5| and having turbulence iB"21|. Un- 



* |galais@ipno.in2p3.fr| 

t ''jpknelle@ncsu.edu' 
1 1 volpe@ipno.in2p3.fr] 



derstanding how the neutrino flavour conversion in su- 
pernovae is modified by the neutrino-neutrino interaction 
is a key theoretical and phenomenological issue. This 
question is at present being intensively investigated. It 
is now clear that three flavour conversion regimes exist. 
They consist first in a collective synchronization of the 
neutrinos with no flavor conversion (2^ [2^ , then in the 
occurrence of " bipolar" flavour oscillations 0, [25| and 
finally in a complete swap of the (anti) neutrino spectra 
above a critical energy (the spectral split) [13, HI, Ha] (for 
a review see [H HI]). 

Several works have investigated the conditions for the 
neutrino-neutrino effects to be triggered and for the on- 
set of the bipolar oscillations. Studies are available both 
in simplified models [l^ and in more sophisticated three- 
fiavour and multi-angle simulations j27l Issj . Using the 
analogy with a pendulum, in p5| it has been pointed 
out that, contrarily e.g. to the well known MSW effect, 
collective fiavour conversion effects occur for any value 
of the third neutrino mixing angle, while such a param- 
eter needs to be stricly non-zero. Besides the authors 
have shown that bipolar oscillations start in inverted neu- 
trino mass hierarchy, because the flavor polarization vec- 
tor is in an unstable position. However, it has recently 
emerged, that such instabilities are a more general be- 
haviour and can be present for any neutrino mass hier- 
archy, depending on the primary neutrino fluxes js^-js^ . 
They can induce not only single but also multiple spec- 
tral splits. Another step in identifying the conditions 
for the bipolar oscillations to start has been done in [23| 
where an heuristic condition involving the vacuum and 
the neutrino-neutrino interaction terms is idcntifled. 

Most of the available calculations of neutrino prop- 
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agation in supernovae including the neutrino-neutrino 
interaction assume the so-called "single-angle" approx- 
imation where the flavour evolution history is trajectory 
independent; while "multi-angle" calculations consider 
the flavour history along different trajectories and dif- 
ferent interaction angles arc considered. While the for- 
mer has been shown to catch well qualitatively and quan- 
titavively many features of the multi-angle calculations 
[soj . the two calculations reveal differences that can be 
important. In particular in multi- angle calculations col- 
lective effects present flavour decoherence, as discussed 
e.g. in (3l| : while, when the matter density exceeds (is 
comparable) the neutrino density, collective effects can 
be strongly suppressed (or be affected by multi-angle 
decoherence) [32[ . In Ref. [231 the authors have also shown 
that the location where bipolar oscillations start varies if 
a single-angle or a multi-angle calculation is performed. 
The implications of different onset locations of the bipo- 
lar oscillations on the r-proccss has been further investi- 
gated in |42| where it was shown that multi-angle versus 
single-angle calculations can produce different r-proccss 
abundances. Note that the importance of the neutrino- 
neutrino interaction for the nucleosynthesis of heavy el- 
ements was identified in an early work [4]| . In the long 
run it is clear that accurate theoretical predictions will 
be based both upon full multi-angle calculations and also 
consider a non-spherical geometry for the neutrinosphcre, 
which can be important in astrophysical environments 
such as accretion disks around black holes and coalescing 
neutron stars [43| . 

In addition to the sensitivity to the details of the neu- 
trino spectra at the neutrinosphcre and the neutrino hi- 
erarchy, the results in presence of the neutrino-neutrino 
interaction can also be sensitive to the Dirac CP violating 
phase, due to loop corrections or to physics beyond the 
Standard Model, as demonstrated in [361, [13 ■ While first 
studied in [s^, the existence and conditions for possible 
CP violatingeffects in supernovae has been established 



in [36[ and [37[ , in presence of the neutrino- neutrino in- 



teraction, and a first quantification of the CP violating 
Dirac phase impact on the neutrino fluxes has been per- 
formed. The analytical results in [36| are now indepen- 
dently confirmed in [s^ . Note that Ref. [1^ has explored 
the impact of a non-zero Dirac phase on the neutrino 
degeneracy parameter at the Big- B ang nucleosynthesis 
epoch, extending the work done in j36l ISTj . In conclu- 
sion, at present, many features of the neutrino-neutrino 
interaction effects have been understood and some of the 
phenomenological implications explored, but many ques- 
tions remain. In particular, work is still needed to fully 
unravel the physical mechanisms underlying the neutrino 
flavor conversion in presence of such contributions, and 
also their interplay with the unknown neutrino proper- 
ties. 

The goal of the present work is to try to gather further 
insight in the neutrino-neutrino effects using the basis of 
the instantaneous eigenstates in matter, instead of the 
flavour one. To this aim we follow the neutrino evolution 



through matter and present general expressions for the 
diagonal and off-diagonal terms of the Hamiltonian de- 
scribing neutrino evolution in the basis of the eigenstates 
in matter, for an arbitrary number of neutrino families 
and for the case of an arbitrary Hamiltonian - in the 
flavor basis - having in particular complex off-diagonal 
contributions. We define the corresponding generalized 
non-adiabaticity parameters. This part extends a previ- 
ous work [3^. While the relations we find are general, 
we apply our results to the case of neutrino evolution 
in core-collapse supernovae, and give explicit expressions 
in the two neutrino flavor case. We then point out for 
the first time the important role of the matter mixing 
matrix Dirac phase, engendered by the presence of the 
neutrino-neutrino interaction contribution to the Hamil- 
tonian. Numerical calculations are provided of the di- 
agonal and off-diagonal terms of the Hamiltonian in the 
basis of the eigenstates in matter. We show that the 
start of the bipolar oscillations is associated with a rapid 
growth of the phase derivative and identify an analytical 
condition that involves the neutrino-neutrino interaction 
and vacuum terms. 

Our paper is organized as follows. The formalism and 
our generalized expressions for the non-adiabaticity pa- 
rameter in the basis of the eigenstates in matter is pre- 
sented in Section II, followed by a focus upon the two 
flavor case as an example. In Section HI we apply our 
results to the case of neutrino propagation in supernova 
where the Hamiltonian is composed of the coupling to 
matter and to neutrinos and present numerical results 
in two flavors for the diagonal and off-diagonal entries 
of the Hamiltonian in the basis of eigenstates in pres- 
ence of matter. Section IV is a conclusion. Appendix A 
shows that the Majorana phases in matter do not change 
our conclusions. Appendix B provides the expressions 
for the derivative of the neutrino-neutrino Hamiltonian 
in the multiangle case. 



II. THE FORMALISM 



In this section we derive general expressions for the 
Hamiltonian, in the basis of the instantaneous eigenstates 
in matter, describing neutrino propagation in an environ- 
ment when the Hamiltonian in the flavour basis contains 
complex potentials. We will call from now on such a ba- 
sis the 'matter' basis. We do not make any assumption 
about the Hamiltonian entries, but only that the Hamil- 
tonian and its derivative can be computed. To make ex- 
pressions more explicit we consider neutrino propagation 
in the framework of core-collapse supernovae. 
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A. Neutrino Evolution in the Flavor Basis 

The neutrino evolution is determined by the Schrodinger 
equation (h = c = 1) : 



dt ^ 



(1) 



where -0'-^^ arc the neutrino amphtudes for a neutrino 
to be in a given flavour state, for an arbitrary number N 
of neutrino families, i?*^-^^ the Hamiltonian in the flavour 
basis. The latter is composed of multiple terms. 



Hif) = UKU'' + H. 



(/) 



Hif) 



(2) 



namely the rotated vacuum Hamiltonian U KU\ the 
i f) 

potentials i?^at due to the coupling of neutrinos with 
matter, which are diagonal in the flavor basis, and the 
neutrino-neutrino interaction term Hi)u which is not, in 
general, diagonal in the flavor basis. The Hamiltonian K 
is given hy K = diag(ki,k2,k^---,kN) with ki being the 
neutrino energy eigenvalues in the mass basis. 

In Eq.([2|) U is the unitary matrix relating the flavor 
and the mass basis, i.e. \va) = X]i=i at t^ail^^i)- For 

= 3 the U matrix is the well known Maki-Nakagawa- 
Sakata-Pontecorvo matrix Q|- Obviously the elements 
of U are related to each other via the requirement of uni- 
tarity. This allows us to express some elements in terms 
of others after specifying the phase of the determinant, 
and restrict the magnitude of the remaining independent 
elements. Thus the N x N mixing matrix U can be 
parametrized in terms of N{N — l)/2 mixing angles and 
N{N +l)/2 phases, out of which TV are of Majorana while 
TV - 1 are of Dirac typeE]. This leaves {N - 1){N - 2)/2 
leftover, CP phases. The dependence upon the Majo- 
rana phases, which we label as a,; can be factored in 
the mixing matrix: i.e. we write U as C/ = A with 

independent of the a's and the matrix A defined to 
he A — diag{cxp{—iai),cxp{~ia2),.-.,cxp{—iaN)) [3- 
The N Majorana phases and the matrix A have no role 
in neutrino oscillations because all observables are 
given by the squared modulus of a matrix element and 
the Majorana phases only ever enter through the phases 
of the elements. 

It is well known that the Dirac phases can be absorbed 
by redefining the charged fermion fields in the standard 
model Lagrangian. This possibility indicates that their 
absolute values cannot affect observables, however, once 
this is done, we have removed this degree of freedom. 
While this rephasing can certainly be performed when 
neutrinos evolve in vacuum, in matter this procedure 
might not be approriate, in particular when the neutrino 
propagation Hamiltonian comprises off-diagonal complex 



terms, as is the case in presence of the contribution com- 
ing from the neutrino-neutrino interaction. Indeed, we 
will show that such phases can play an important role. 

If we restrict to iV = 3 families and consider that neu- 
trinos interact via the standard weak interaction with 
ordinary matter composed of electrons, protons and neu- 
trons, the contribution to H^^'> coming from H^l^ is di- 
agonal in the flavor basis y*^-^-* (r) = ^^^(^^(r), 0, 0) with 
Ve(r) = \/2GFne{r) and ne(r) the electron density, r 
being the distance in the supernova. The contribution 
coming from scattering on neutrons can be substracted. 
In fact it is possible to remove the trace of this Hamilto- 
nian because we are free to add or subtract an arbitrary 
multiple of the unit matrix (including a term that is a 
function of position) , since the only effect of such a term 
is to introduce an overall phase. The third term in H^^^ 
correspond to the neutrino-neutrino interaction Hamilto- 
nian j23| 

H^fJ = V2Gf f - q • ^)dn,Jq' (3) 

where dn^^ [dnp^ ) are the differential neutrino (anti- 
neutrino) number densities and p^^ {Pva ) is the density 
matrix for neutrinos (anti-neutrinos) , enconding neutrino 
flavour conversion, whose expression is e.g. for neutrinos 



\Ve? Ve^l VeK 
12 



(4) 



and similalrly for . In Eq. Q the neutrino- neutrino 
interaction term is built from the contributions coming 
from neutrinos and anti-neutrinos having different mo- 
menta q, q', q = q/|ql with and born with the flavour a 
at the neutrinosphero The fact that their flavour might 
have changed up to the interaction point is encoded in 
the neutrino density matrix. 



B. Neutrino Evolution in the 'Matter' Basis 

We now consider the basis that instantcnously diagonal- 
izes the Hamiltonian H^^'> Eq.® 



c/t ii = k 



/fci 

fc2 



(5) 



where fc^, i = 1, are the corresponding energy eigenval- 
ues in matter. For the two and three flavor cases, the 



^ The assignment of the phases to Dirac or Majorana type is not 
unique: different parameterizations will change the assignment. 



^ The fact that the neutrino flavour content might have changed 
is indicated with a flavour underlined. 
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expressions for the eigenvalues are well known [46|, |47[ 
and one may also find an algebraic expression for TV = 4. 
For N > 5 the Abel-Ruffini theorem indicates that no 
general algebraic formula exists. Once the eigenvalues 
are found, one can determine the unitary transformation 
U relating the flavour basis and the 'matter' basis. Note 
that, to distinguish quantities evaluated in vacuum from 
those calculated in matter, in the following we will use 
symbols with tilde, to signify that e.g. the energy eigen- 
values as well as the parameters (phases and angles) of 
the U matrix are those evaluated in presence of matter. 
Like U, the matrix U also has TV Major ana phases cti 
that we can factor out from U so that U = ttfA with 
A = diag{exp{—iai) , exp(— za2), , exp(— laAr)). If the 
eigenvalues are non-degenerate then we can derive the 
identity that the elements of U satisfy: 



c. 



Ill 



(i) 

a/3 



(6) 



fact, if we differentiate the eigen-equation H'^f)u = ilk 
and multiply the result by J7^, we derive the result that 



dt 



The commutator vanishes for the diagonal elements of 
this equation, because K is diagonal. Thus we find the 
result that 



a,l3 



dt 



ill) 



By considering the off-diagonal elements for this same 
expression, the following equation is deduced 



dt 



dt 



0, i^] (12) 



where C'*-' is the co-factor matrix of i?^-^' — ki and so 



c 



/3/3 



Z^7 ^77 



(7) 



Note that both equations are independent of the a's. The 
denominator of Eqs.(IB][7]) is the trace of the cof actor ma- 
trix H^^^ — ki which is basis independent, so that 

7 j 

In practice we use Eq.® to relate the elements in a col- 
umn to one specific element f/^i, then use Eq.([71) to eval- 
uate c7^J3 

Once U is found, a change of basis can be made using 
[/, to give the Schrodinger equation in the 'matter' basis 



j— - 

dt 



K 



dt } ^ 



(9a) 
(9b) 



where H is the Hamiltonian in the basis of the eigenstates 
in matter. The term [/^ dU/dt appears because the mat- 
ter eigenvalues are functions of time, which requires that 
U also be a function of time. In order to evaluate H 
we need to compute dU jdt. The derivation of the 
off-diagonal elements of U'^ dU j dt is straight-forward. In 



Again, since K is diagonal, the off-diagonal elements of 
the commutator have a compact form : 




1 

6ki 



where Skij is Skij ^ ki — kj . 



dt 



(13) 



The remaining missing pieces of H are the diagonal 
terms of dU /dt. Because of the Majorana phase ambi- 
guity mentioned previously there is no unique expression 
for these terms but, by using Eq.®, we can eliminate the 
ambiguity from all elements of a column of U bar one, 
Upi. Using the identity that C/t dU/dt = -dUydtU we 
derive that 




u: 



c4{ 



Upi Uq 



-iu*. 

"'dt 



/3i 



u, 



dU, 



Pi' 



dt 



dU, 



dt 



(14) 



The first term of Eq. p^ is well defined because it is inde- 
pendent of the a's; while the Majorana phase ambiguity 
is entirely contained in the second term. Introducing the 
variable Qi to represent the first term in Eg. p^ . and 
using Eq.®, we find 



iQ^ =^2^ 



c. 



dt 



c 



13a 



dt 



(15) 



^ Note that since the Majorana phases on do not appear in 
Eqs.| |6I7I |. we have no way to determine them and can make any 
choice including functions of the position (see Appendix A). 



withS = (2|f//3^|2 (Tr(C«))^)-i. This leaves the second 
term in Ea. ((T4| . If the phase of Ufa is chosen to be solely 
the Majorana phase i.e argiJJgi) = —ai - since there are 
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TV — 1 Dirac phases /3 will be fixed to the electron flavour 
for all i - then 



Inserting this equation into H we obtain our final result 
that it can be written as 



1 



2\Ufi,\ 



Pi 



01' 



dt 



dt 



duj 



Thus we obtain the result that the diagonal elements of 
i/t dU/dt are 




dai 
~dt 



(17) 



H = 



ki +Qi 



' dt 



: {5ki2+SQi2) 



2 TT 



^ 12 
13 



{Ski2+SQl2) 



ri 



k2 + Q2-^^ 



dt 

■ {Sk23+6Q23) p 

27r ^ 



23 



■ (3fci3+i5Qi3) p 

' 27r ^13 

■ {Sk2i+SQ2-6) p 

* 2^ ^23 

^■3 + Q3-^^ 



(18) 



with 5Qij given by 5Qij = Qi ^ Qj and we have intro- 
duced the Tij functions which are the generalized non- 
adiabaticity parameters for neutrino oscillations with ar- 
bitrary potentials. These quantities are defined as 



2 7re**"'^ 
{6kij + 6Qij 

2 Tre"'"'^ 



Skij {Skij + 5Qij) 




(19) 



W — - — U 



where we have substituted in, the results from Eg. ([131). 
The non-adiabaticity paramaters Fy have, in general, 
both real and imaginary components whose origins are 
from different components in the flavor-basis Hamilto- 
nian. Let us investigate this connection further by con- 
sidering the case of two-neutrino flavors. While many 
aspects are familiar from the MSW problem, new compo- 
nents arise when we consider more general Hamiltonians. 



C. The application to two-neutrino flavors 

Let us consider the general parametrization of the mixing 
matrix in matter for two neutrino flavours 



U 



1 




cos f sm U 
— sin 9 cos 9 







(20) 



and note that the phase /3 is the Dirac phase in matter. 
With this parametrization the elements of H^^^^ are given 



(21) 



ki cos^ 9 + k2 sin^ 
-Ski2 e 



ip 



cosf smf 



-5ki2 e *^ cos 6* sm9 
ki sin^ 9 + k2 cos^ 9 



From Eqs. (|2T|) it is very simple to derive several dif- 
ferent relationships between the matter angle 9 and the 
elements of the flavor basis Hamiltonian e.g. 



sm 



Hee - fci 

k2 - ki 



(22) 



The expression for cos^ 9 is very similar to equation ((22)) 
but with Hfj_f^ in place of Hee- For the Dirac phase we 
derive that 



tan/3 



3(geM) 
^H(Fe^) 



(23) 



From this equation we immediately see that in the case 
when the off-diagonal elements of H'^f^ are independent 
of position then (3 is simply the same as its vacuum value. 
By differentiating U and multiplying by f/'l' one find^ 



dU_ 
'dt 



sin 20 — i5qi2 



sin 20 gti55i2 



cos 




e 

-,-?(5qi2 







d9 
'dt 



(24) 



* Note how the Majorana phases do not appear in these expres- 
sions. 

^ We have omitted the derivatives of the Majorana phases. 
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so that, following Eq. p^ . the full Hamiltonian 



H 



V 



ki +Qi 



2tt 



12 




(25) 



explicitly read^ 



H = 



fci + /3 sin^ 



iSai2 sin 29 



+ it 



128 



(26) 

where 5ai2 ~ ai — 012- From the diagonal elements of H 
we read off the Q's to be 



Qi =/3sin2 
Q2= P cos^ 



(27a) 
(27b) 



and from the off-diagonal entries wc derive that the gen- 
eralised non-adiabaticity parameter is 



ri2 = --F 



27re* 



sin 26* ~;\ 



',12 



and observe that it depends both on the derivative of 
the matter angle 6, as in the MSW case, and on the 
derivative of the matter phase From Ea.ipS]) one sees 
that the are defined up to the Majorana phases. In 
order for the imaginary component of T12 to be non-zero 
we require that the off-diagonal elements of TJ^-^^ rotate 
in the Argand plane. Neutrino-neutrino interaction gives 
exactly such a term because the flavor basis Hamiltonian 
includes contributions from the density matrices i.e. H ^ 
P- 

We can then differentiate Eas. (PT|) and eventually find 



= cos^ 6 Hee + sin^ 6 



-cos 61 sine* (e^^iJ, 



\.e) (29a) 



k2 ~ SIT? 9 Hee + COS^ 6 Hn 



cose sin^ ( e^'^Hea + e~''^i/„ 



(29b) 



which is consistent with our expectations from Ea.(|lip. 
We can then differentiate the expressions for and (3 and 
derive that 



iifee - i^MM ) ^fcl2 - {Hee - H^^) Ski2 



(30) 



and 



sm20iHee-H^ 
25ki2 



cos 261 



25k 



12 



(31) 



From now on wc indicate differentiation by / = df/dt for com- 
pactness. 



The first term in Eq. pip will be familiar to many readers 
because it is the MSW term, the second term is new. 

From the two terms on the r.h.s. of Eq. (|23p . two ex- 
pressions for /3 can be obtained : 



Ski2 sin 20 



(32a) 



(32b) 



If one substitutes into Eq.([28|) our results for 6 and /3, 
one gets 



l^cos 6 sin 6 {jle. 



ee H i_Li_L 



Ski2{dki2 + SQ12) 

+Het, cos^ e e"^ - H^e sin^ 9 e^^'^) (33) 



The reader may verify that this is exactly the expression 
one would have obtained directly from equation ([Tl 



III. THE APPLICATION TO THE 
NEUTRINO-NEUTRINO INTERACTION 
PROBLEM 

A. The phase derivative and the onset of the 
synchronization and bipolar regimes 

Let us now apply the formalism and the results of the 
previous section to the Hamiltonian Eqs.([2]|3]) describing 
neutrino propagation in a core-collapse supernova with 
the neutrino-neutrino interaction contribution. In the 
calculations we use a realistic matter density profile com- 
ing from supernova simulations. We treat the neutrino- 
neutrino contribution in the single-angle approximation. 
In Appendix B we provide equations for the derivative 
of the neutrino-neutrino interaction Hamiltonian in the 
multi-angle case. Note that from now on we replace time 
(t) with distance (x). If we calculate the derivative of 
H necessary to determine the diagonal and off-diagonal 
terms of the Hamiltonian in the 'matter' basis one has 
two contributions, namely H = Hmat + Hpy^ since the 
vacuum term does not contribute. 

In the single- angle approximation, H^i, Eq.Q reads 



Tjs.a. 



with the geometrical factor 



F(x)^"- 



y/2GF 1 
27ri?2 2 



1 



Ry 

X 



(34) 



(35) 



where is the radius of the neutrinosphere. The non- 
linear contribution is given by 

a •' 

(36) 
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with the neutrino flux at the same location for a 

neutrino of a flavor a. The corresponding derivative Hin^ 
includes contributions from both the derivative of the 
geometrical factor and the density matrices, i.e. 

Hl-;^- ^ Fixy-^-Gipy-"- + F(a;)"-"-G(p)"-"- (37) 

where the calculation of F{x) is straightforward and the 
one for the non-linear term is obtained using of the Li- 
ouville Von-Neumann equation. This gives 

G(p)-- = ~iJ2 [ {[H,p,Jq')]L,Jq') 

+ [H,p^Jq')]*L^^iq'))dq' (38) 

with H indicating the Hamiltonian for anti-neutrinos. In 
the case of anti-neutrinos, Eg. ([55)1 holds but by replacing 
H ^ H and ^ i>. 

Using the results just derived one can identify two in- 
teresting conditions to be fullfilled : i) for the neutrino- 
neutrino interaction effects to occur; ii) for the onset of 
bipolar oscillations. Let us first focus on the the start of 

the synchronization regime and discuss the value of /? at 
X = Riy. In particular, for such a quantity to be non zero, 
it is necessary Tjg^ 7^ and not parallel to He/j, ■ If H^fj, is 
pure real then i?e^ has to have an imaginary component. 
Using Eqs. p4j|35p and ^ one obtains 



ac 



H, Pv^ 



TTvac 



(39) 
(40) 



and the same for anti-neutrinos but with H instead of 
H. By taking the explicit vacuum terms 



TTVac 



TTvac 

"f" ~ 4q' 



sin 29 



one gets 



J dq' ^sin20 J dq'-^ 



(41) 



(42) 



If we replace the neutrino number flux at the neutri- 
nosphere by e.g. the Fermi-Dirac distribution for the 
j/g contribution, one obtains 



I dq' [H,p,^]^^L,^ 



Am^Lo . (77)^1(77) 
■ sm 20- 



(43) 

with Fi{ri) the complete Fermi-Dirac integrals (without 
the gamma function normalization) and 



1 



1 



(44) 



where {Ey) indicates the neutrino average energies and 



47ri?2 



-Lo- 



(45) 



It has been pointed out in Ref. [25| , in the analogy with a 
pendulum, and also in Ref. j30| that the vacuum mixing 
angle need to be non-zero for the neutrino-neutrino in- 
teraction to have an effect. In our case, where we exactly 
solve the two-neutrino evolution with all contributions 
in the 'matter' basis, since the dominating term comes 

from /3 we find that the condition to have effects com- 
ing from the neutrino-neutrino interaction needs Eas. (|42p 
and (|33]) to be non-zero, i.e. if sin 26* 7^ 0, in agreement 
with Ref. m, 113. 

Let us now focus on the onset of bipolar oscillations. In 
Ref. [131, within a three flavor multiangle calculation, it 
has been pointed out heuristically that the transition to 
the bipolar regime should be due to a condition involving 
the vacuum and the neutrino-neutrino interaction contri- 
butions. Here with only two flavors a condition signifying 
the end of the synchronization regime can be explicitly 

identified from j3 Ea. p2bp namely that |i?e;iP has to ap- 
proach zero. More explicitly, the element _ffe^t only in- 
volves the vacuum and self-interaction terms so we find 
the onset of the bipolar regime occurs when: 



1 2 I Tjvac I 2 
- \J^eii I " 



l-^e^ I 



2 n{H:^''H*n ^ (46) 



This is agreement with the heuristic condition given in 
[27| . Note that, although obtained in two flavours, our 
condition pS)) is obtained from the equations derived in 
section II. C where no simplifying assumption has been 
made regarding the neutrino Hamiltonian. 



B. Numerical results and discussion 

Here we present our two-neutrino flavor numerical results 
in the 'matter' basis obtained using Eas. (P51) . (|32b[) . ([571 - 

In our calculations the parameters are fixed as fol- 
lows. For the mixing parameters we take \dm'^\ = 
2.4 X 10~^eV^, 6* = 9°. For the matter density profile 
we use pb = 1-5 x lO'^ (10/a;) in units of g.cm""^ and of 
km (for x). The neutrinosphere Ri, is taken at 10 km and 
the corresponding neutrino fluxes are assumed to be of 
Fermi-Dirac type with average energies {E^J = 12 MeV, 
{Ei,J = 15 MeV and {E^J = 18 MeV. Equipartition of 
energy is assumed, with a total luminosity of 4 x 10^^ 
erg.s"^, so that in our calculations collective effects only 
appear in inverted hierarchy. However our findings con- 
cerning the role of the matter phase are expected to hold 
also in the case where equipartition is not satisfied and 
multiple spectral splits appear. 

Figure [1] shows the oscillation probabilities of the neu- 
trino instantenous eigenstates in matter, as a function of 
the distance in the supernova for two different neutrino 
energies. The case considered here is inverted hierar- 
chy. The synchronization regime, bipolar oscillations and 
the spectral split are easily recognized. The first corre- 
sponds to the region where neutrino fiavour conversion is 
frozen; while in the second regime oscillations appear in 
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FIG. 1. Two-flavor neutrino evolution in the 'matter' basis: 
The flgures show the neutrino (left) and anti-neutrino (right 
flgures) survival probabilities for first (solid black) and second 
(solid grey) matter eigenstates as a function of distance within 
a core-collapse supernova for energies of 5 MeV (upper) and 
10 MeV (lower) figures. 
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FIG. 3. Same as Figure [2] but for the off-diagonal element 
of two-neutrino flavor Hamiltonian in the basis of the in- 
stantenous eigenstates in matter. The results correspond 
to the real (upper), imaginary (middle) parts of and 
its modulus |ri2| (lower flgures), multiplied by the factor 
{Ski2 + 5Q 12) /5~ki2 ■ 



Neutrinos Anti-neutrinos 




20 40 60 80 100 120 140 160 180 20 40 60 80 100 120 140 160 180 200 



r [km] ;■ [km] 

FIG. 2. Two-flavors neutrino evolution in the 'matter' ba- 
sis: Diagonal elements of the Hamiltonian Eq. (|26|) (fci + 
Qi)/5ki2 (upper), {k2 + Q2)/5ki2 (middle) and their differ- 
ence 1 + 5Q\2/5k\2 (lower flgures). The curves correspond to 
a 5 (black) and 10 MeV (grey) neutrino energy for neutrinos 
(left) and anti- neutrinos (right flgures). The black (grey) lines 
show the average for a 5 (10) MeV neutrino. The calculations 
include the vacuum mixing, the coupling to matter and the 
neutrino-neutrino interaction. 



the probabilities. The spectral split phenomenon occurs 
between 70 and 100 km and correspond to the complete 
flavour conversion, giving rise to a swap of the neutrino 
spectra, depending on the neutrino energies. As an ex- 
ample, the two energies shown are smaller (larger) than 
the spectral critical energy for which no (full) neutrino 
flavour conversion takes place. 



Figure [2] presents the diagonal [ki + Qij /5ki2 with 
i = 1,2 calculating using Egs. iP^ and (|32b|) . as well as 
the difference of the diagonal elements 1 + 5Qi2/6ki2. It 
is indeed this quantity which is important to follow the 
neutrino flavor evolution. The results arc shown both 
for neutrinos and anti-neutrinos and different neutrino 
energies. One can see that, since 9 is minimum (w 0) for 
neutrinos and maximum (« ti'/Z) for anti- neutrinos, one 
of the diagonal matrix elements is approximately given 
by hi (ki for neutrinos and ^2 for anti-neutrinos) while 

the other oscillates very fast due to /?. 

The real and imaginary part of the off-diagonal contri- 
butions Ti2{Ski2+SQi2)/6ki2 Eq. ipS]) of the Hamiltonian 
in the 'matter' basis Eg. ((26)) are presented in Figure [S] 
Results are given for both neutrinos and anti-neutrinos 
of different energies. Comparing the scales of Figures [5] 
andOone sees that the diagonal elements are much larger 
than the off-diagonal ones as expected. The off-diagonal 
contributions are practically zero in the synchronization 
region below 50 km where they abruptly become non- 
zero. Fast oscillations coming from the derivative of the 
phase are again present. Note that the numerical results 
are very similar both for neutrinos and anti-neutrinos. 
Figure [3] presents the ratio of the average of the modulus 
of the numerator over the average of the denominator of 
the generalized adiabaticity parameter Ea. ([28|) . One 
can see that such a quantity is larger than one in the re- 
gion where the neutrino-neutrino interaction effects are 
dominant. 

Let us now discuss the relative role of the diagonal and 
off-diagonal terms of the Hamiltonian in the 'matter' ba- 
sis, when the neutrino-neutrino contribution is included. 
One first notes that, in absence of such a contribution, 
SQi2 = and Ski2 would not vary over the region of 200 
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km near the neutrinosphere. However, after the inclusion 
of H^^, the difference between the diagonal entries in H, 
i.e. Ski2 + SQ 12, drops suddenly at around 50 km when 
the bipolar oscillations begin (bottom part of Figure [2]) . 
This is due to the fact that one of the two diagonal el- 
ements suddenly changes due to the P cos^ 9 (/? sin^ 6) 
for neutrinos (anti-neutrinos). As a consequence, when 
bipolar oscillations start, the difference between the di- 
agonal elements oscillates around zero, due to the differ- 
ence SQi2 = /3 cos 261. These results make it clear that 
the Dirac phase /3 acquires a particularly significant role 
in the neutrino flavor evolution. It leads to an oscillatory 
degeneracy between the diagonal elements of H, increas- 
ing significantly the importance of the off-diagonal entries 
in H which had a negligible role during the synchroniza- 
tion regime. Indeed, in the region where the neutrino- 
neutrino interaction effects are dominant the derivative 
of the matter mixing angle is 6* « 0. The size of the off- 
diagonal elements of H is approximately (59. It is clear 
that these small off-diagonal contributions become very 
important when the difference of the diagonal elements 
falls to zero. 



Neutrinos Anti-neutrinos 
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^47 48 49 50 51 52 47 48 49 50 51 52 



/■ [km] r [km] 

FIG. 5. Contributions to the off-diagonal matrix elements 
|7/ep| Eq. (|46p of the two- neutrino flavor Hamiltonian. The 
curves show : (dashed), jHe/l'l^ (black) and the total 

l-ffeuP (grey), in units of eV^ , for neutrinos (left) and anti- 
neutrinos (right). Results for a 5 IVIeV neutrino are given in 
the upper figures and for a 10 MeV one in the lower figures. 



40 60 80 100 120 40 60 80 100 120 140 

r [km] r \kin] 

FIG. 4. Ratio of the average of the moduhis of the numer- 
ator over the average of the denominator of the generalized 
adiabaticity parameter Eg. psp corresponding to the two- 
neutrino flavor Hamiltonian in the 'matter' basis, for neutri- 
nos (left) and for anti- neutrinos (right figure), with an energy 
of 5 (black) and 10 IVleV (grey). 



Concerning the onset of bipolar oscillations. Figure [5] 
shows that the different terms corresponding to the con- 
dition |i?e/iP = Eq. (|46| . At the beginning iJg^ is dom- 
inated by the constant vacuum contribution. Since the 
modulus of l^fg^P increases as the interference between 
these two terms. Only when these contributions have 
been cancelled off can the phase (3 vary rapidly. So the 
point where \H^I^'^\ = \H^'^\ engenders an abrupt change 
in Qi and Q2 (Figure [2]) and leads to a concomitant sud- 
den increase in ri2. When the neutrino- neutrino inter- 
action term becomes negligible another change in the Qi 
occurs. From Figure[2]this appears to occur at an energy 
dependent distance of 110 km for 5 MeV and 130 km for 
10 MeV neutrinos, where the neutrino probabilities have 
reached their asymptotic behavior. 



IV. CONCLUSIONS 



General expressions for the diagonal and off-diagonal 
terms of the Hamiltonian in the basis of the instantenous 
cigenstatcs in matter have been derived, corresponding to 
a general Hamiltonian in the flavour basis with complex 
contributions. Relations have been given to determine 
in particular the off-diagonal contributions as a function 
of the neutrino mixing matrix relating the flavor basis 
and the one formed by instantenous eigenstates in mat- 
ter. We have applied these findings to the case of neu- 
trino propagating in a core-collapse supernova in the two- 
ncutrino case. Our analytical results show that the flavor 
evolution is governed by the derivatives of both the mix- 
ing angle and a phase. The numerical results show the 
important role of the matter phase. In particular, the 
diagonal and off-diagonal contributions to the Hamilto- 
nian, in the 'matter'basis, show abrupt changes at the 
end of the synchronization region that we have associ- 
ated to a divergence condition for the derivative of the 
matter phase, involving the vacuum and the matter term. 
Such a condition, that characterizes the onset of bipolar 
oscillations, is only slightly dependent on the neutrino 
energy. In future work we plan to extend the present 
study to the case of three flavors, in order to determine 
whether additional insight may be gathered. 
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V. APPENDIX A 



VI. APPENDIX B 



Here wc show that the Majorana phases in matter do not 
influence the neutrino flavour conversion. In fact, there 
exists a basis in which the contributions to the diagonal 
matrix elements of the matter Hamiltonian are absent. 
This occurs if we introduce the adiabatic basis, discussed 
in [39! , related to the matter basis via the unitary trans- 
formation W{x) i.e. ■0 = W ^p^'^K In this new basis, the 
Schrodinger equation is 



dt 



dt 



dt 



(47a) 
(47b) 



The matrix W is chosen so that it removes the diag- 
onal elements of including the derivatives of the 
Majorana phases that appear in the diagonal entries of 
dU /dt. This requirement indicates that W is diagonal 
so writing W as 



/ exp[— 2i7r0i] 



W = 







exp[— 2 i7r(/)2] 



(48) 



V 



one obtains that the phases 0i are defined by 



^ _ J_ f I ,^ doi 
dt " 27r r' + ^' dt 



Integrating gives 



1 

27 



-&^{t) + / dt' (h + Qt 



(49) 



(50) 



Thus wc find that W factors into a matrix ]^ inde- 
pendent of the a's. and the Majorana matrix A, i.e. 
W = \^ — A'^ . When we insert this solution for W 
into H^""^ we find that the Majorana phase dependence 
in the off-diagonal entries of disappears too i.e. 



dt 



dt 



(51) 



Thus we are able to find a basis in which the Hamiltonian 
is entirely independent of the Majorana phases which, in 
turn, implies the S-matrix 5'^°' is also independent of 
the phases in this basis. The difference between matter 
basis and this new basis is simply a position dependent 
rephasing of the elements in the Hamiltonians and the 
associated S matrices. Thus we reach the expected con- 
clusion that that the Majorana phases we introduced in 
the matter mixing matrix cannot affect neutrino oscilla- 
tion observables regardless of what we pick for them. As 
a consequence, in our derivations, we can set the Majo- 
rana phases to constants (or zero) and neglect them (as 
well as their derivatives), without any loss of generality. 



We furnish here expressions for the derivative of the 
neutrino- neutrino interaction Hamiltonian Eq.Q, that 
can be of use in future multi-angle calculations of neu- 
trino evolution in a core-collapse supernova, using, as 
in the present manuscript, the 'matter' basis. Assuming 
spherical geometry for the neutrinosphere, the non-linear 
Hamiltonian Eq.Q reads 



V2Gf 



TTTn.a. 

" 2TTm 



E 



JcOSe^ax 

[p.M',o')L.M')~plM.o')L,^{q') 



dq' dcos9' (1 — cos 6* cos 6*') 
(52) 



with: 



cos 9r, 



'l-{R.lxY 
Using the Leibniz Integral Rule : 



(53) 



- f{y,z)dt = 

CZ Ja(z) Ja{z) 



dz 



dt 



+f{b{z).z) 
-f{a{z),z) 



dbjz) 

dz 
da{z) 

dz 



(54) 



The derivative for Eq. (f52|) is 
^/2GF 



Tjra.i 

with 

R{e\q') 



2nm 



dq' {R{e',q')-S{x)T{q')) 



(55) 



S{x) 



dcosO' {—i (1 — cos 6* cos 6*') 
+ [H,p,Jq',0')]* L,Jq')] 



(56a) 
(56b) 



x^l~{R,/xf 

T{q') = (1 - COs6'cOs6'„iaa;) [PvAq' ■,Omax)L^J^q') 

-Pp„('?'i^'mQx)ip„ ((?')] (56c) 
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